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Abstract 


We  develop  an  approximation  and  convergence  theory  /or  Galerkin  approximations 
to  infinite  dimensional  operator  Riccati  differential  equations  formulated  in  the  space  of 
Hilbert- Schmidt  operators  on  a  separable  Hilbert  space.  We  trea/  the  Riccati  equation  as 
a  nonlinear  evolution  equation  with  dynamics  described  by  a  nonlinear  monotone  pertur¬ 
bation  of  a  strongly  coercive  linear  operator.  We-prcn^e  a  generic  approximation  result/or 
quasi-autonomous  nonlinear  evolution  systems  involving  accretive  operators  which  wg  then 
usedto  demonstrate  the  Hilbert-Schmidt  norm  convergence  of  Galerkin  approximations  to 
the  solution  of  the  Riccati  equation.  •We-rHtrstrate  the  application  of  our  resultsjn  the 
context  of  a  linear  quadratic  optimal  control  problem  for  a  one  dimensional  heat  equation. 


*This  research  supported  in  part  by  the  United  States  Air  force  Office  of  Scientific 
Research  under  grant  AFOSR-87-0356.  Part  of  this  research  was  carried  out  while  the 
author  was  a  visiting  scientist  and  consultant  at  the  Institute  for  Computer  Applications 
in  Science  and  Engineering  at  the  NASA  Langley  Research  Center.  Hampton,  VA.  which 
is  operated  undei  NASA  contract  NASlTSiU?. 
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1.  Introduction 

Recently  a  great  deal  of  attention  has  been  given  to  the  analysis  of,  and  the  develop¬ 
ment  of  approximation  theory  for,  infinite  dimensional  operator  Riccati  differential  equa¬ 
tions.  This  is  due  in  part  to  the  role  they  play  in  the  characterization  of  the  closed-loop 
feedback  structure  of  the  solution  to  linear-quadratic  optimal  control  problems  for  dis¬ 
tributed  parameter  systems  (see  [10]).  Working  from  Temam’s  [14],  [15]  formulation  (see 
also  Barbu  [3])  of  a  class  of  operator  Riccati  equations  in  the  space  of  Hilbert-Schmidt  op¬ 
erators  on  a  separable  Hilbert  space,  we  develop  an  approximation  and  convergence  theory 
for  generic  Galerkin  type  approximations.  We  make  the  usual  assumption  that  our  un¬ 
derlying  Hilbert  space  admits  a  densely,  continuously,  and  compactly  embedded  subspace. 
The  only  condition  that  we  then  require  on  our  sequence  of  approximating  Galerkin  sub¬ 
spaces  is  the  usual  one  that  the  corresponding  sequence  of  orthogonal  projections  converge 
strongly  to  the  identity  in  the  stronger  topology.  We  are  able  to  obtain  Hilbert-Schmidt 
operator  norm  convergence  of  the  approximating  Riccati  operators,  uniformly  in  time  on 
compact  time  intervals. 

Our  effort  here  differs  significantly  from  other  recent  treatments  of  the  approximation 
theory  for  infinite  dimensional  operator  Riccati  equations  (see,  for  example,  [G])  in  that,  i) 
we  obtain  Hilbert-Schmidt  as  opposed  to  strong  operator  convergence,  and,  more  impor¬ 
tantly,  ii)  our  convergence  theory  is  based  directly  upon,  and  involves  only,  the  differential 
equation  itself  rather  than  equivalent  integral  equations.  Of  course  in  order  to  do  this 
we  must  necessarily  consider  a  somewhat  more  restrictive,  but  still  sufficiently  interesting 
from  an  applications  point  of  view,  class  of  problems.  For  example  we  require  that  the 
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linear  part  of  the  equation  be  strongly  coercive  on  a  space  of  Hilbert-Schmidt  operators 
and  that  the  non-homogeneous  or  quasi  autonomous  perturbation  (the  state  penalization 
operator  in  the  cdntext  of  the  LQ  control  problem)  be  Hilbert-Schmidt.  Our  treatment 
here  is  restricted  to  the  quasi-autonomous  or  constant  coefficient  case.  The  temporally 
inhomogeneous  problem  requires  a  different  approach;  our  results  for  the  time  dependent 
case  will  be  reported  elsewhere. 

Our  convergence  theory  is  based  upon  a  generic  approximation  result  for  nonlinear 
quasi-autonomous  evolution  equations  in  Banach  space  with  dynamics  described  by  accre¬ 
tive  operators.  In  section  2  we  prove  a  nonlinear  analog  of  the  well-known  Trotter-Kato 
theorem  on  the  approximation  of  linear  semigroups  -  i.e.  stability  (uniform  dissipativity) 
and  consistency  (strong  operator  convergence  of  the  resolvent)  yield  convergence  (see,  for 
example,  [11]).  The  result  that  we  prove  here  is  closely  related  to  similar  approximation 
results  for  nonlinear  evolution  systems  which  have  appeared  elsewhere  in  the  literature  (for 
example,  [4]  and  [8]).  However  the  theorem  we  prove  below  is  one  that  is  most  appropriate 
for  the  particular  class  of  problems  that  are  of  interest  to  us  here.  We  follow  Goldstein 
[8]  and  give  a  proof  using  an  idea  originally  suggested  by  Kisynski  [9]  for  the  linear  case 
wherein  convergence  is  demonstrated  via  an  application  of  an  existence  theorem  to  an 
evolution  equation  in  an  appropriately  constructed  sequence  space. 

In  section  3  we  briefly  outline  Temam's  [15]  formualtion  of  the  Riccati  equation  as  a 
well  posed  nonlinear  evolution  equation  in  the  space  of  Hilbert-Schmidt  operators.  Our 
approximation  and  convergence  theory  is  developed  and  presented  in  section  4.  In  section 
5  we  illustrate  the  application  of  our  results  in  the  context  of  a  linear  quadratic  optimal 
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control  problem  (i.e.  the  linear  quadratic  regulator  problem)  for  a  one  dimensional  heat 
equation. 


2.  An  Abstract  Approximation  Result  for  Nonlinear  Evolution  Equations 

Let  Xo  be  a  real  Banach  space  with  norm  denoted  by  j  •  j0.  Let  XJ  be  its  dual  and 

let  Aq  :  Xo  — ►  2Xo  be  an,  in  general,  multi-valued  nonlinear  closed  accretive  operator  (i.e. 

{(x,  y)  :  x  e  X0,  y  e  A0x}  is  a  closed  subset  of  Xo  xX0,  and  |x]  — x2|0  <  |xi  +  Ay!  —  (x2  + Ay2)|0 

for  all  A  >  0,  Xj  £  Xo,  and  yj  e  AoXj,  i  =  1,2)  on  Xo-  Define  the  domain  of  A0  to  be  the  set 

Dom(Ao)  =  {x  £  Xo  :  Aox  ^  <p},  and  the  range  of  Ao  to  be  the  set  7?.(Ao)  =  (J  Aox. 

xeDom(Ao) 

Since  A0  is  accretive,  it  follows  that  for  A  >  0  the  resolvent  of  Ao  at  A,  J0(  A )  =  (I  +  AAo)-1, 
is  a  well  defined  single-valued  nonexpansive  operator  (i.e.|J0( A )y i  —  Jo(A)y2|o  <  |y j  —  y2|o, 
yi,y2  e  Dom(J0(A)))  defined  on  Dom(J0(A))  =  7£(I  +  AA0).  Suppose  that  T  >  0,  let 
t  — >  fo(t)  be  an  X0  —  valued  map  defined  on  [0,T]  and  let  x£  e  X0.  We  consider  the  initial 
value  problem  in  X0  given  by 

(2.1)  x0(t)  +  A0x0(t)  B  f0(t).  a.e.  t  z  (0, T), 

(2.2)  xO(0)  =  Xq. 

We  shall  say  that  a  function  xo(-)  :  [0.T]  — ►  Xo  is  a  strong  solution  to  the  initial  value 
problem  (2.1),  (2.2)  if  it  is  continuous  on  [0,T],  Lipschitz  on  every  compact  subinterval  of 
(0,T),  differentiable  almost  everywhere  on  (0.T),  and  satisfies  (2.1)  and  (2.2).  We  shall 
call  xo  an  integral  solution  of  (2.1 ).  (2.2)  it  is  continuous  on  [0.T1  satisfies  r2.2),  and  if 
the  inequality 

<  fo(r)  ~  y ,x0(r)  -  x  >0  dr 


-|x0(t)-x|jj  <  ~ jx0 ( s ) 


X  0  + 


/ 
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holds  for  all  x  e  Dom(Ao),  y  e  Aox,  and  0  <  s  <  t  <  T  where  the  pairing  <  •  >o-' 

X0  x  Xo  — *  R  is  given  by  <  y,x  >0=  sup{x*(y),x*  e  3>0(x)}  with  <£0  :  X0  — ►  2X°  denoting 
the  duality  map  on  Xo  (see  [3]).  A  strong  solution  of  (2.1 ),  (2.2)  is  of  course  also  an  integral 
solution.  The  following  existence,  uniqueness,  and  regularity  results  are  given  by  Barbu 
[3], 

Theorem  2.1  Let  {X,  |  •  j }  be  a  real  Banach  space  and  suppose  that  A  :  X  — ►  2X  is  a 
closed  accretive  operator  on  X.  Let  C  be  a  closed  convex  cone  in  X  such  that  Dom(A)  C  C 


and  7Z(l  +  AA)  D  C,  A  >  0.  If  x°  e  Dom(  A)  and  f  e  Lj(0,  T;  X)  with  f(t)  e  C,  a.e.  t  e  (0,  T), 
then  there  exists  a  unique  integral  solution  x  to  the  initial  value  problem 


(2.3) 


(2.4) 


x(t)  +  Ax(t)  5  f(t),  a.e.  t  e  (0,T), 


x(0)  =  x°, 


with  x(t)  e  Dom(A)  a.e.  t  e  (0,T).  Moreover,  if  y  is  the  integral  solution  to  (2.3),  (2.4) 
with  f  replaced  by  g  e  Li(0,T;X),  then 

|x(t)  -  y(t)|2  <  |x(s)  -  y(s)|2  +  2  [  <  f (r)  -  g(r),x(r)  -  y(r)  >  dr 


for  0  <  s  <  t  <  T. 


Theorem  2.2  If  in  addition  to  the  hypotheses  and  conditions  of  the  previous  theorem  we 
have  x°  e  Dom(A)  and  f  e  T:  X).  then  the  initial  value  problem  (2.3),  (2.4)  has  a 

unique  strong  solution  x  with  x  e  W1  ^(0,  T;  X).  x(t)  e  Dom(A)  a.e.  t  e  (0,T)  and 

|x(t)|  =  |Ax(t)  +  f(t)|  <  |Ax°  +  f(0)|  +2  f  |f'(s)jds,  a.e.  te(0,T). 

Jo 
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Furthermore,  if  X  and  X*  are  uniformly  convex,  then  x  is  everywhere  differentiable  from 
the  right  and 

£±(t)  =  (Ax(t)  +  f(t))°,  0  <t<T 

where  (Ax(t)  +  f(t))°  denotes  the  unique  element  of  minimum  norm  in  the  set  Ax(t)  -fi 
f(t). 

Henceforth,  we  shall  assume  that  there  exists  a  closed  convex  cone  Co  C  Xo  for  which 
Dom(Ao)  C  Co  and  TZ(1  +  AAo)  D  Co,  A  >0,  that  Xq  e  Dom(  A0 )  and  that  f0  e  Li(0,  T;  X0 ) 
with  fo(t)  e  Co  for  a.e.  t  e  (0,  T).  Consequently  Theorem  2.1  applies  and  we  are  guaranteed 
that  the  initial  value  problem  in  Xo,  (2.1),  (2.2),  admits  a  unique  integral  solution. 

We  prove  an  approximation  result  in  the  spirit  of  the  well  known  Trotter-Kato  theorem 
for  the  approximation  of  linear  semigroups  (see,  for  example,  [11]),  and  the  approximation 
theorems  given  by  Crandall  and  Pazy  [4]  and  Goldstein  [S]  in  a  nonlinear  setting.  Although 
the  convergence  theorem  we  shall  prove  here  does  not  differ  significantly  from  the  ones  given 
in  [4]  and  [8],  the  latter  results  are  stated  inappropriately  and  are  somewhat  too  restrictive 
for  the  application  we  intend  to  consider  below.  Following  Goldstein  [8],  the  proof  we  shall 
give  here  is  based  upon  an  idea  first  suggested  by  Kisynski  [9]  in  the  context  of  linear 
semigroup  approximaton.  We  argue  convergence  via  an  application  of  the  existence  and 
uniqueness  result.  Theorem  2.1. 

For  each  n  =  1,2,...  let  Xn  be  a  closed  linear  subspace  of  Xo.  and  let  A„  :  Xn  —*  2Xn  be 
a  closed  accretive  operator  on  X„.  As  was  the  case  above  with  Aq.  for  A  >  0  the  resolvent 
of  A„,  J„(A)  =  (I  4-  AA„)_I  is  a  well  defined,  single-valued,  nonexpansive  operator  defined 
on  Dom(J„(A))  =  7Z(I  -f  AA1().  We  assume  that  there  exists  a  closed  convex  cone  Cn  in  Xn 


for  which  Dom(A„)  C  C„  and  71(1  +  AA„ )  D  C„,  A  >  0.  We  assume  further  that  for  each 
n  =  1,2,  is  an  element  in  Dom(A„)  and  that  f„  is  a  function  in  Lj(0,  T;X„)  with 

fn(t)  e  Cn  for  almost  every  t  e  (0,T).  It  then  follows  that  Theorem  2.1  implies  that  the 
initial  value  problem  in  Xn  given  by 

(2.5)  x„(t)  +  A„xn(t)  3  f„(t),  a.e.  t  e  (0,T), 

(2.6)  x„(0)  =  x„, 

admits  a  unique  integral  solution  xu  e  C(0,T;X„). 

In  our  discussions  below,  we  shall  use  the  notation  lim  Dn  D  Do  where  Dn  C  Xn  and 

n — ►  oo 

Do  C  Xo-  By  this  we  shall  mean  that  for  each  zq  e  Do  there  exists  a  sequence  {zn}^j  with 

z„  f  D„  and  lim  zn  =  z0.  Our  fundamental  approximation  and  convergence  result  is  given 
n— • ‘oo 

in  the  following  theorem. 

Theorem  2.3  Suppose 

(i)  lim  Cn  D  C’o 

n  — ♦oo 

(ii)  lim  fn(t)  =  fo(t)  and  there  exists  age  Li(0,T)  for  which  |fn(t)|  <  g(t),  n  =  1,2,..., 

n  — -oo 

for  almost  every  t  e  (0,T) 

(iii)  limn-ocX°  =  x° 

(iv)  lim  Jn(A)yn  =  Jo(^)yo  for  each  A  >  0  whenever  y„  e  Cn  and  limn_00yn  =  y0  f  C0. 

n— *oo 

Then  lim  xn  —  x0  in  C(0,T;Xo)  where  x„  for  n  =  1,2,...  and  x0  are  the  unique  integral 

n  — *oo 

solutions  to  the  initial  value  problems  (2.5),  (2.6)  and  (2.1),  (2.2)  respectively. 

Proof:  If  we  define  the  linear  space  X  over  the  reals  by 

X  =  {u  =  {u„}“0  :  u„  e  Xn,n  =  0, 1,  2, ...,  lim  un  =  u0} 

n  —  oo 
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and  for  u  =  {un}^0  e  X  set  |u|  =  sup  |un|o,  then  {X,  |  ■  |}  is  a  real  Banach  space.  Let  C 

II 

be  the  closed  convex  cone  in  X  given  by  C  =  {u  =  {u„ },^0  e  X  :  un  e  C„,n  =  0,1,2..} 
Note  that  if  Co  ^  <f>,  hypothesis  (i)  implies  that  C  ^  <?.  Define  the  operator  A  :  X  — ►  2X  by 


Au  = 


—  {^n}n=0  ‘  ^ 

n  =  o,  l,2,...,lim„_00  vn  =  e0},  u  =  {«„}“0  e  Dom(A) 
4>  u  £  Dom(A) 


Dom(A)  =  {u  =  {un  }~0  e  X  :  u„  e  Dom(  A„ ),  n  =  0, 1, 2  . . . ,  and  for  each  n  =  1,2,... 
there  exists  a  v„  e  Anun  for  which  lim  vn  =  v0  e  A0u0}. 

n  —  cc 

The  operators  An  being  closed  and  accretive  implies  that  the  operator  A  is  closed  and 
accretive  as  well.  Indeed,  for  example,  for  A  >  0,  u*  e  Dom(A),  i  =  1,2,  and  v‘  e  Au1 ,  i  =  1, 2, 
we  have 


lu1  -  u2|  =  sup|ui  -  u2|0  <  sup  |uj  +  Av*  -  (u2  +  Av2)|0  =  |uJ  +  Av1  -  (u2  +  Av2)|. 

n  n 

The  fact  that  A  is  dosed  can  be  argued  analogously. 

Clearly  Dom(A)  C  C  and  it  also  follows  that  for  A  >  0  1Z(1  +  A  A)  e  C.  To  see  this  let 


A  >  0  and  let  v  =  {vn}^0  e  C.  Then  vn  e  C„  and  Cn  C  7?.(I  +  AAn),  n  =  0, 1,2, ...  imply 
that  we  may  define  un  =  J„(A)vn,  n  =  0,1,2,..  and  set  u  =  {u„}£T0.  Now  veCcX  and 
hypothesis  (iv)  yield  lim  u„  =  lim  J„(A)vn  =  J0(A)v0  =  u0;  thus  u  e  X.  Also,  if  we  set 

n — oc  n— *oo 

Wn  =  (vn  — un)/A,  n  =  0,1,2,..  then  w  =  {wn}^L0  e  X,  w„  e  A„u„  and  lim  wn  =  w0  e  A0u0. 

n  — -oo 

It  follows  that  un  e  Dom(  A„ )  and  therefore  that  u  e  Dom(  A)  with  w  e  Au  and  v  e  (1  + AA)u. 


We  conclude  that  v  e  7J(I  +  AA)  and  therefore  that  C  C  7?.(I  -f  AA). 

Define  z°  =  {x°}^i0.  Then  hypothesis  (iii)  implies  that  z  e  X.  Moreover,  using  (iv) 
it  can  be  shown  that  Dom(A)  =  {u  =  {un}£L0  e  X  :  u„  e  Dom(An),n  =  0,1,2,...}.  Thus 
z  6  Dom(A).  If  we  define  the  function  f  by  f(t)  =  {f„(t)}^0  for  almost  every  t  e  (0,  T),  then 
hypothesis  (ii)  implies  that  f  e  Li(0.  T;X),  and  f„(t)  e  C’„,  n  =  0.1,2,...  a.e.  t.  f  (0.T)  implies 
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that  f(t)  e  C,  a.e.  t  e  (0, T).  It  follows  from  the  auguments  given  above  and  Theorem  2.1 
that  the  initial  value  problem  in  X 

z  +  Az(t)  3  f(t),  a.e.  t  e  (0,T) 

z(0)  =  z° 

admits  a  unique  integral  solution  z  e  C(0,T;X)  with  z(t)  e  Dom(A)  for  almost  every 
t  e  (0,  T).  In  addition,  inspection  of  the  proof  of  Theorem  III.  2.1  given  in  [3]  (paying 
particular  attention  to  how  solutions  are  actually  constructed)  reveals  that  we  must  have 
z(t)  =  {x„(t)}£T0  where  for  n  =  1,2,...,  x„  is  the  unique  integral  solution  to  the  ini¬ 
tial  value  problem  (2.5),  (2.6)  in  Xn  and  xo  is  the  unique  integral  solution  to  the  initial 
value  problem  (2.1),  (2.2)  in  Xo-  Consequently,  since  z(t)  e  X,  t  e  [0, T],  it  follows  that 
lim  x„(t)  =  x0(t)  for  each  t  e  [0,T], 

n — ►oo 

Finally,  that  the  convergence  of  x„(t)  to  x0(t)  is  in  fact  uniform  in  t  for  t  e  [0,T]  can 
be  argued  as  it  was  done  by  Goldstein  in  (S|.  Let  e  >  0  and  t  e  [0,T]  be  given.  Then 
z  =  {x„}^L0  e  C(0,  T;X)  implies  that  |z(t)  —  z(s)|  <  e/3  for  all  0  <  s  <  T  with  |t  —  s|  < 
26  =  26(e).  Now  lim  xn(s)  =  xo(s)  for  each  s  e  [0,T]  implies  that  |xn(s)  —  xo(s)|o  <  e/3 

n — *-00 

for  all  n  >  Ns  =  N(s,e).  Thus 

|xn(t)  -  x0(t)|0  <  jxn( t)  -  xn(s)|0  +  |x„(s)  -  x0(s)|0  +|x0(s)  -  x0(t)|0  <  e, 

for  all  0  <  s  <  T  with  j t  —  s|  <  26  and  n  >  Ns.  Let  k  be  the  greatest  integer  less  than  or 
equal  to  T/6  and  set  N  =  max{N0,  N«,  N2<i,  •  • . ,  N^^,  Nt|.  Then  |xn(t)  —  x0(t)|o  <  e  for  all 
n  >  N  =  N(e)  and  the  proof  of  the  theorem  is  complete. 
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3.  Operator  Riccati  Equations  on  Spaces  of  Hilbert-Sehmidt  Operators 

In  this  section  we  briefly  review  and  outline  Temam’s  [15]  results  on  operator  Riccati 
equations  set  in  the  Hilbert  space  of  Schmidt  class  operators  on  a  separable  Hilbert  space. 
Let  H  be  a  real  separable  Hilbert  space  with  inner  product  and  associated  induced  norm 
denoted  by  (•,  •)  and  |  ■  |  respectively.  Let  V  be  another  real  separable  Hilbert  space 
with  inner  product  <  •,•  >  and  induced  norm  j|  •  |j  .  We  assume  that  V  is  densely  and 
continuously  embedded  in  H  (i.e.  V  =  H,  and  there  exists  a  positive  constant  p  for  which 
M  <  p||<p!|,<p  e  V).  Identifying  H  with  its  dual.  H*,  we  have  V  t— ►  H  =  H*  <— »  V*  with 
the  embedding  of  H*  in  V*  being  dense  and  continuous  as  well.  Denote  the  usual  operator 
norm  on  V*  by  ||  ■  ||».  Let  7  denote  the  canonical  isomorphism  (Riesz  map)  from  V  onto 
V*.  Then  for  <p,  xj>  e  V  we  have 

(777,  V’)  =<  VM/’  > 

where  (•,•)  in  the  above  expression  denotes  the  usual  extension  of  the  H  inner  product  to 
the  duality  pairing  between  V*  and  V.  We  assume  further  that  the  embedding  of  V  into 
H  is  compact.  It  then  follows  that  7-1  e  £(V,  V*)  D  £(H,  V)  and  that  j~1  is  self-adjoint, 
positive,  and  compact  as  a  mapping  from  H  into  H.  We  note  that  V*  is  in  fact  a  Hilbert 
space  with  inner  product  <•,•>»  given  by 

<  v?,i/’  >.  =  <  7~V^~V  >  =  (y>,7-V),  (  v*. 

We  have  that  ||(p||,  =  s/<  ^  f  V'. 

Since  7-1  is  self-adjoint,  positive,  and  compact  on  H  there  exists  an  orthonormal  basis, 

{ek }ic=  1 1  f°r  H  such  that  7-1ek  =  k  =  1,2,...  for  some  real  numbers  p^,  k  =  1, 2 . 

Consequently  {pj^ekjj^-j  and  {pkCklkii  are  orthonormal  bases  for  V  and  V*  respectively. 


9 


For  separable  Hilbert  spaces  X  and  Y,  let  HS(X,Y)  denote  the  Hilbert  space  of  Hilbert- 
Schmidt  operators  from  X  into  Y.  Denote  the  corresponding  inner  product  and  induced 
norm  by  [-,-]hs(X,Y)  ancl  1  '  Ihs(X.Y)  respectively.  Set  H  -  HS(H,H)  with 

'X) 

[$,*]*  =  [^,^]ns(H,H)  =  5j$ek,tfek). 

k=  1 

and  !$!•«  =  |$|hs(H,H)  =  €  H.  Define  the  Hilbert  space  V  by  V  = 

HS(V*,H)  fl  HS(H,  V)  with  inner  product 

OO  OO 

[$,*]v  =  [$i  ^]hs(V*,H)  +  [$,  ^]hs(ii,V)  =  ^2pli^k,9ek)  +  ^2  <  ^ek,^ek  >, 

k=  1  k=l 

|j <E> ||  v  =  $]y,  e  V.  It  is  not  difficult  to  show  that  the  dense  and  continuous 

inclusions  HS(V*,H)  C  HS(H,H)  C  HS(V,H),  and  HS(H,V)  C  HS(H,H)  C  HS(H,V*) 
hold.  Also,  it  can  be  argued  that  HS(V*,H)  and  HS(V,H),  and  HS(H,V)  and  HS(H,V*) 
are  dual  pairs  with  respect  to  the  duality  pairing 

OO 

[$,'!']*  =  5>ck,*ck). 
k=  1 

It  follows  that  V*  =  HS(V,H)  +  HS(H,V*)  and.  identifying  H  with  its  dual  H* ,  that 
V  <— ♦  H  <— >  V*,  with  the  inclusions  dense  and  continuous. 

Let  a(-,-)  :  V  x  V  — >  R  be  a  bounded  strongly  V-elliptic  bilinear  form  on  V  x 
V.  That  is,  there  exist  constants  a.  ,3  >  0  for  which  >  o||^>||2  and  |a(^,t/')l  < 

/^IMIII^II*  Vi  4’  e  V.  (With  minor  modification  our  theory  continues  to  hold  if  the  form 
a(-,  )  satisfies  the  coercive  inequality  a(y?,9>)  +  A|^|2  >  a||^||2,  \p  e  V,  for  some  A  f  R. 
To  keep  our  presentation  here  as  simple  as  possible,  however,  in  our  discussions  below  we 
shall  treat  only  the  case  A  =  0.)  Let  A  e  £(V,Y*)  be  the  operator  defined  by  (Ad) A')  — 
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a (0,0),  <t>,4’  e  V.  Let  a*(-,-)  :  V  x  V  — >  R  be  the  form  which  is  adjoint  to  a(-,-)<’  that  is 
a*(<i,V’)  =  a(V’. <?),  f  V.  Then  a*(.p,^)  >  o Jj^Jj2  and  |a*(^,^)l  <  /?IMII|£||,  <P,V  c  V. 
Let  A  e  £(V,V*)  he  the  operator  defined  by  =  a *(r>,v)  =  a(i/\y>),  e  V. 

Let  C0  be  the  closed  convex  cone  in  H  given  by  C’o  =  {*F  e  HS(H,H)  :  $  =  >  0}. 

Let  T  >  0  and  suppose  that  llo  e  Co  and  that  Fo  e  Li(0,T;H)  with  Fo(t)  e  Co  for  almost 
every  t  e  (0,T)  are  given.  Let  <F  — ►  Jro(<5)  be  a  single  valued  map  defined  for  each  $  f  Co 
with  range  in  Ti  which  is  continuous  from  7i  into  itself.  Assume  further  that  To  has  the 
property  that 

[^0(*)-.Fo(n*-*]w  >  o, 

for  all  $,4/  (  Co-  We  note  that  if  z  — >  JF0(z)  is  a  single  valued  complex  function  of  the 
complex  variable  z  with  J’o(O)  =  0  and  which  is  analytic  on  the  nonnegative  real  axis,  then 
the  mapping  $  — +  satisfies  Tq( $ )  e  H  for  all  $  e  Cq  and  is  continuous  from  'H.  into 

'H  (see  Dunford  and  Schwartz  [5],  Theorem  XI. 6. 7. 7).  We  seek  a  solution  ITo  to  the  initial 
value  problem 

(3.1)  rio(f)  +  >cn0(t)  +  n0(tM  +  ^o(no(t))  =  Fo(t),  a.e.  t  (  (0,  T) 

(3.2)  n0(0)  =  n« 

with  ITo ( t )  (  Co.  a.e.  t  e  (0,  T). 

Note  that  when  To($)  —  <F2,  (3.1)  becomes  the  standard  quadratic  Riccati  equation. 
Elementary  properties  of  Hilbert-Schmidt  operators  (see  [14])  can  be  used  to  argue  that 
for  <F,  *F  e  Co  we  have 

[<F2  -  T2,  <F  -  <F]W  =  [<F(<F  -  \F).  <F  -  \F]«  +  [(<F  -  'F)'F.<F  -  $]«  >  0. 


Essentially  the  same  argument  can  be  used  for  the  case  of  /o(<F)  =  <Fn  where  n  is  any 


positive  integer.  For  more  general  quadratic  terms,  for  example,  the  one  that  results  in 
the  case  of  the  linear-quadratic  optimal  control  problem,  ^"o(^)  =  with  M  e  £(H), 

M  =  M*,  and  M  >  0,  once  again  a  similar  argument  will  work  so  long  as  >  0.  $  e  Co 
(see  Proposition  2.2  in  [15]).  This  will  of  course  be  true  if  for  all  $  e  Co 

Remark  3.1  At  this  point  a  comment  regarding  the  relationship  between  the  operator 
A *  and  the  adjoint  of  A  is  in  order.  Define  the  operator  A  :  Dom(A)  C  H  — ►  H  to 
be  the  restriction  of  the  operator  A  to  the  set  Dom(A)  =  {y?  e  V  :  Ap>  e  H}.  It  can 
be  shown  (see  [12])  that  Dom(A)  =  H  and  consequently  that  A  admits  an  H-adjoint, 
A*  :  Dom(A*)  C  H  — *  H.  The  operator  A*  is  the  extension  of  A*  to  an  operator  defined 
on  all  of  V  or,  equivalently,  A*  is  the  restriction  of  A *  to  the  H-dense  subset  Dom(A*)  = 
{<p  t  V  :  A*p>  e  H}.  We  note  also  that  -A  is  the  infinitesimal  generator  of  a  uniformly 
exponentially  stable  semigroup,  {T(t)  :  t  >  0},of  bounded  linear  operators  on  H.  Similarly, 
—  A*  is  the  infinitesimal  generator  of  the  adjoint  semigroup  |T(t)*  :  t  >  0}  on  H.  In 
addition,  it  can  be  argued  that  both  of  these  semigroups  admit  respectively  restrictions 
and  extensions  to  analytic  semigroups  on  V  and  V*  (see  [1],  [13]). 

An  appropriate  refomulation  of  (3.1 ),  (3.2)  will  allow  an  application  Theorems  2.1  and 
2.2.  Define  the  operator  Lo  (  C(  V.  V'  )  by  Lo^  =  +  4>«4.  for  <J>  e  V.  It  is  not  difficult 

to  argue  that  Lq  is  strongly  V  —  elliptic:  that  is,  there  exists  a  constant  uj  >  0  for  which 

(3.3)  [L0$.$U  >c’||$||2v, 

and  therefore  that  the  set  Dom(Lo)  ={4*6  V:  Lo'h  f  H }  is  dense  in  7i  (see  [13]).  Define 


the  operator  Ao  :  Dom(  Ao  )  C  %  — ♦  “H  by 

A0$  =  L0$  +  ^o(^),  $  e  Dom(A0)  =  Dom(L0)  n  C0. 

It  follows  that  Ao  is  a  closed  operator  on  7i  and  that  it  is  strongly  V  —  monotone.  That  is, 
for  f  Dom( Ao)  we  have 

[Ao<&-Ao'J',<&-,l']*  =  [L0($-#),$-$]h  +  [^o($)-.^o(’I,),$-’I']-h  > 

From  this  it  can  be  argued  at  once  that  Ao  is  accretive  on  its  domain  and  that  7^(1  + AAq)  D 
Co  for  all  A  >  0. 

We  rewrite  (3.1),  (3.2)  as  the  initial  value  problem  in  H  given  by 

(3.4)  iio(t)  +  A0Ilo(t)  =  F0(t),  a.e.  t  e  (0, T) 

(3.5)  n0(0)  =  n° 

Now  Dom(A0)  =  Dom(Lo)nCo  =  C0  and  recall  that  it  was  assumed  that  Dio  e  Co-  Also 
it  was  assumed  that  F0  e  L1(0,T;7i)  with  Fo(t)  e  Co,  a.e.  t  e  (0,T).  Consequently  Theorem 
2.1  yields  the  existence  of  a  unique  integral  solution  Do  to  the  initial  value  problem  (3.1), 
(3.2)  with  n0(t)  e  Co  for  almost  every  t  e  (0,  T).  If  it  is  further  assumed  that  e  Dom(Ao) 
(i.e.  that  ng  e  C0  and  that  .4‘llg  +  e  HS(H,H))  and  that  F0  e  H)  with 

Fo(t)  e  Co,  a.e.  t  e  (0,T),  then  Theorem  2.2  implies  the  existence  of  unique  strong  solution 
Do  e  W1,oo(0, T;  H)  with  A*Do(t)  +  Do(t)A  e  HS(H,  H),  a.e.  t  e  (0,  T).  Using  density  it  can 
be  argued  further  that  if  e  Co  and  Fo  e  L2(0,T;  V*)  with  Fo(t)*  =  F0(t),  F0(t)  >  0,  a.e. 
t  e  (0,  T),  then  there  exists  a  unique  solution  n0  e  L2(0,  T;  V)(l  C( 0.  T;  H)  D  H^O,  T;  V*) 
with  n0(t)  e  Co  for  every  t  e  [0,  T]  (of  course  in  this  setting.  Do  is  a  solution  of  (3.4),  (3.5) 
in  a  distributional  or  V*  sense). 
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4.  Approximation  Theory 

For  each  n  =  1,2,...  let  Hn  be  a  finite  dimensional  subspace  of  H  with  Hn  C  V  for  all 
n.  Let  Pn  :  H  — ►  Hn  denote  the  orthogonal  projection  of  H  onto  Hn  with  respect  to  the 
standard  inner  product,  (•,•),  on  H.  We  shall  require  the  assumption  that 

(4.1)  lim  1 1 P n<p  —  y?||  =  0,  f  e  V. 

n— *oo 

Note  that  (4.1)  implies  that  lim  |Pn>?  —  ~p\  —  0,  ^  e  H  and  that  Pn  is  uniformly 

n  — oo 

bounded  in  n  in  both  £(H)  mid  £(V)  with  respect  to  the  uniform  operator  topologies. 

It  will  be  necessary  for  us  to  define  an  extension  of  the  operator  P„  to  a  bounded 
operator  from  V*  into  V*.  We  do  this  as  follows.  For  <^>  e  V*  define  P„ <^>  e  Hn  to  be  the 
representer  in  Hn  of  the  functional  i p  restricted  to  a  functional  on  Hn.  That  is  P„v?  =  <pn 
where  ipn  is  the  unique  element  in  Hn  (guaranteed  to  exist  by  the  Riesz  representation 
theorem  applied  to  the  Hilbert  space  {Hn,  (•,•)})  which  satisfies  (p,9n)  —  (ipn  ,0n),  0neHn. 
Since  Hn  C  V  C  V*,  Pn  can  be  considered  to  be  a  linear  map  of  V*  into  itself.  Indeed 
PnV7  =  ‘r’n  is  the  continuous  linear  functional  on  V  given  by  (P,,^.#)  =  (^>n,  9), 9  e  V. 
The  definition  of  the  orthogonal  projection  of  y?  e  H  onto  Hn,  (Pn^  -  'r’A)  =  0Jnf  H n , 
reveals  that  we  have  in  fact  defined  an  extension  of  P„  to  an  operator  on  V*  with  range 
Hn  considered  as  a  subspace  of  V*.  Using  the  fact  that  for  ^  e  V*  and  v  e  V  we  have 
( P  n  r5 ,  V’ )  =  (Pn'r’iV’)  ~  (Pn^-  V’  -  Pnf)  =  (Pn>^,  Pnv)  =  (^,  PnU),  and  that the  Pn  are 
uniformly  bounded  in  £(V).  it  is  not  difficult  to  show  that  the  P„  are  uniformly  bounded  in 
£(V*).  It  is  worth  noting  that  the  extension  of  P„  defined  above  agrees  with  the  operator 
that  would  be  obtained  if  we  were  to  extend  Pn  in  the  usual  way  by  considering  it  as  a 
bounded  operator  defined  on  the  dense  subset  H  of  the  Banach  space  V*. 
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It  follows  from  assumption  (4.1)  that  lim  ||Pn^  —  'r’ll*  =  0,  <p  e  V*.  Indeed,  for 

n  — *oc 

^  e  H  we  have  jj Pn v3  —  ^11*  ^  MlPnv*  ~  v5!  ~ *  0  as  n  — *  oo.  Since  H  is  dense  in  V*,  and 
the  Pn  are  uniformly  bounded  in  £(V*),  it  follows  that  ||Pn<^  —  'r’ll*  - 1 *  0  as  n  — >  oo  for 
<p  e  V*.  For  ip,  ip  e  V*,  we  have 

<  Pnp,\l>  >,  =  (Pn ¥Mr-V)  =  (<r\Pn7~V)  =  (^,7_SPn7_V)  =  <  V?i7Pn7~V  >.  • 

Thus  P*  e  £(V*),  the  adjoint  of  P„  considered  as  a  bounded  operator  on  V*,  is  given  by 
P*  =  7Pn7-1-  Assumption  (4.1)  yields 

IIP;*  -  *9||.  =  l|7~1P:*-7_V||  =  ||7"17P»7_,*-7"l*ll  =  ||Pn7_V  -  7"V||  —  0 

as  n  — >  oo,  for  each  <p  e  V* . 

For  each  n  =  1,2,...  define  the  finite  dimensional  subspace  7fn  of  H  by 

Hn  =  {$nPn  :  e  £(Hn)}. 

Note  that  Hn  is  clearly  a  subspace  of  both  H  and  V  since  Hn  finite  dimensional  implies 
that  all  operators  in  Hn  are  of  finite  rank.  Define  the  closed  convex  cone  Cn  C  Hn  by 
Cn  =  {4>„Pn  c  Hn  :  4>n  =  $*,$„  >  0}.  Using  the  fact  that  Hn  C  V  and  Cn  C  C0,  we 
define  the  operator  An  :  Dom(An)  C  Hn  — *  Hn  by 

(4.2)  An($nPn)  =  {A0($nPn)}Un,$nP„  e  Dom(An)  =  Cn. 

That  is  An($nPn)  is  the  element  Ao(4>nPn)  in  V*  restricted  to  a  linear  functional  on  7in. 
Since  Hn  is  a  finite  dimensional  Hilbert  space,  the  Riesz  representation  theorem  implies 
that  An(<£nPn)  =  4' n P „  e  H„  where  'I'aPn  is  the  unique  element  in  Hn  which  satisfies 

[Ao($nPn),  ©nPnIw  =  [^nPni  ©iiPn]^,  ©nPn  C  Hn  ■ 
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The  operators  An  defined  in  (4.2)  are  in  effect  the  standard  Galerkin  approximation  to 
A0.  Furthermore,  the  definition  (4.2)  leads  to  the  same  approximation  to  A0  that  would 
be  obtained  via  the  more  conventional  procedure  wherein  the  operator  A  in  the  definition 
of  Lo  is  replaced  by  its  Galerkin  approximation.  Indeed,  for  each  n  =  1,2,...  define 
the  operator  Au  e  C(Hn)  by  An<p „  =  .Av’niHn,  ft  f  Hn.  That  is,  .An^n  =  V’n  where 
V’n  is  that  element  in  Hn  (once  again  whose  existence  and  uniqueness  is  guaranteed  by 
the  Riesz  representation  theorem)  which  satisfies  (,A^n,#n)  =  (t/’n,#n)i  #n  f  H„.  Noting 
that  since  P„  is  the  orthogonal  projection  of  H  onto  Hn,  it  is  not  difficult  to  argue  that 
[$P„, On]x  =  [$, 0n]x  for  ad  $  e  V*  and  0n  e  Hn ,  it  then  follows  that  for  $nP„  and 
* nPn  6  Hn  we  have 

[An(^nPn)i^uPnlw  —  [Ao(^nPn)t  ^nPn]H  =  [Ao(  A*nPn  )Pn ,  'Pn  Pn]-H 
=  [A**n  P„  +  ^„P„^Pn+/-0($„Pn)P„,^nPn]7i 

oo 

=  52  {(-A*^nPnek,^nPnek)  +  ($ „Pn v4Pnek ,  Vn Pnek )  +  (^b($nPn)Pnek,  'I'nPnek)} 

k=l 

OO 

=  E{M*^nPnek,^„Pnek)  +  (APnek,  ^nPnek )  +  (Pb  ($nP„)Pnek,  ^P^k)} 

k=l 

oo 

=  E{(^$nPnek,^nPnek)  +  (Ai  Pn^k,  A>*^nPnek)  +  (Pn^b($nPn )Pnek ,  'I'nPnek)} 

k=l 

OO 

=  £  ({A*n  +  +  P„^o(^nPn)}Pnek,  ^nPnek); 

*=1 

or 

A„($nPn)  =  {A$n  +  $.,A  +  A(  A>„Pn  )}Pn , 

where  ^(A*)  =  Pn.Po($),  $  e  Co-  In  particular  when  Fo($)  =  $2,  we  have  An($nP„)  = 
{A*n$n  +  <t>nAn  +  *l}Pn. 

For  each  n  =  1,2,...  define  11°  e  Cn  by  11°  =  PnIIgPn  and  Fn  e  Li(0,  T;Wn)  by 
F„(t)  =  P„Fo(t)Pn,  for  almost  every  t  e  (0,T).  We  consider  the  problem  of  finding  a 


solution  n„  to  the  initial  value  problem  in  Hu  given  by 


(4.3)  li„(t)  + Ann„(t)  =  F„(t),  a.e.  t  e  (0,T) 

(4.4)  nn(o)  =  n°n 

with  n„(t)  e  Cn  for  almost  every  t  «  (0,T).  The  definition  of  the  operator  A„  together  with 
the  properties  of  the  operator  Ao  yield  that  for  each  n  =  1,2,...  A„  is  closed  and  strongly 
V  —  monotone  on  its  domain.  Hence  An  is  accretive  with  7^(1  +  AAn)  D  Cn  for  all  A  >  0. 
Since  Dom(An)  =  Cn  =  Cn,n°  e  Cn,  and  Fn  e  Li(0,T;Hn)  with  F„(t)  e  Cn  for  almost 
every  t  e  (0,T),  Theorem  2.1  yields  the  existence  of  a  unque  integral  solution  Iln  to  the 
initial  value  problem  (4.3),  (4.4)  with  n„(t)  t  Cn,  a.e.  t  e  (0,T). 

We  shall  argue  the  convergence  of  nn  to  IIo  as  n  — ►  oo  in  C(0,  T;  H)  (i.e.  the  Hilbert 
Schmidt  norm  convergence  of  n„(t)  to  Ho(t),  uniformly  in  t  for  t  e  [0,  T])  via  an  application 
of  Theorem  2.3.  In  order  to  do  this  we  shall  require  some  preliminary  lemmas.  The  first 
lemma  below  is  a  technical  lemma  which  can  also  be  found  in  [7].  For  completeness  we 
state  it  here  and  have  included  its  rather  brief  proof. 

Lemma  4.1  If  {aj}?^  is  an  absolutely  summable  sequence  of  real  numbers,  then  there 

exist  sequences  {b,}-^.  and  {ci}^,  such  that  lim  bj  =  0,  {cj}^,  is  absolutely  summable, 

i— *oo 

and  a;  =  b;c;,  i  =  1,2,...  . 

OO 

Proof:  Let  a  =  ^  Jaj  |  and  for  j  =  0,1,2,..  define  the  nonnegative  integers  kj  as  follows. 
i  =  l 

Let  ko  =  0  and  let  kj  denote  the  first  index  for  which 
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Set  b;  =  1/j  and  c;  =  ja;,  for  i  =  kj-i  +  j  =  1,2,...  .  Then  b;c;  =  a;, 

CO  CO  OO 

i  =  1,2, ...,  lim  bj  =  0,  and  53  lcil  =  53  J  S3  lak  I  <  a  +  53  tt  <  °o. 

i=l  j  =  l  k=kJ_1  +  l  j  =  l 

Lemma  4.2  Let  X  and  Y  be  real  separable  Hilbert  spaces  with  inner  products  denoted  by 
<  •,  •  >x  and  <  •,  •  >y ,  respectively.  If  $  e  HS(X,  Y)  then  $  can  be  factored  as 
with  e  C( Y)  compact  and  <£2  e  HS(X,Y). 


Proof:  Let  {.r,}^,  and  {y,}^  be  orthonormal  bases  for  X  and  Y  respectively.  Since 
$  e  HS(X,Y)  it  has  a  representation  in  the  form  of  an  infinitie  matrix  $  <-»  =  [< 

OO  OO  OO 

yj,$Xj  >y]  with  53  53  ^1}  <  oo.  For  i  =  1,2,...  set  a;  =  53‘r’lj-  Since  the  sequence 

.=  1.7  =  1  j  =  i 

{ai}i^i  is  absolutely  summable,  we  can  apply  Lemma  4.1  and  obtain  sequences  {bj};Ti 

OO  OO 

and  {ci}|Ti  for  which  aj  =  b,Ci,i  =  1,2, ....  lim  bj  =  0,  and  53  ci  =  53  lc>l  < 

‘~*00  »=i  j^i 

OO 

oo.  Define  e  £(Y)  and  $2  e  £(X,Y)  by  ^y  =53  v^7  <  y,y.  >r  y.,  y  c  Y,  and 

t=i 

OO  OO 

4>2x  =  53  13  <  x’xi  >X  yi,  x  e  X.  Then  $!$2  =  <3>,  and  since  lim  \/bt  —  0 

i=i  j=i 

/  \  2  oo  oo  oo  oo 

and  53  ~  13  F  53  'filj  =  53  —  53  c»  <  °°’  it  follows  that  <3>1  is  compact  and 

i,j  '  '  '  1=1  ’  J  =  1  «=1  *  «=1 


$2  e  HS(X,  Y). 


Lemma  4.3 

(a)  lim  |P„^Pn  -  =0,  $  f  H. 

n  — *oo 

(b)  lim  ||Pn$Pn  —  $j|v  =  0,  $eV. 

n  — ►oo 

Proof:  (a).  We  consider  Pn  to  be  an  element  in  £(H).  Then  for  $  e  7i  —  HS(H,H)  we 

have 

|Pn$P„-$| n  <  |Pn$Pn~Pn$|*  +  |Pn$-*| n 
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<  |<£Pn  ~  $|-W  +  |Pn$  — 

=  |($P  „)*-$*!*  +  |P„*-$|* 

=  |Pn$*  -$*|*  +  |Pn$-$|*. 

Since  $,<£*  e  HS(H,H),  we  can  apply  Lemma  4.2  with  X  =  Y  =  H  to  obtain  $  =  $1$2 
and  $*  =  ($*)1($*)2  with  $1,($*)1  e  £(H)  compact  and  $2,($*)2  e  HS(H,H).  Then 

|P„$  —  =  \{P„  ~  /)$  $  IhS(W,H)  —  ~  I}®  L(H)  l$  I HS(HM)  ' 

The  fact  that  P„  converges  strongly  to  the  identity  on  H  and  that  e  £(H)  is  compact 
imply  |Pn<£  —  d>|"H  — ♦  0  as  n  — 1  oo.  Similar  estimates  for  |P„$*  —  yield  the  desired 

result. 

(b).  Let  $  e  V  =  HS(V*, H)  D  HS(H,V)  and  consider  ||Pn$P„  -  $||HS(V.H)- 
Since  |Pn|  =  1  it  follows  that 

||Pn$Pn  -  $I|hS(V,H)  <  l|Pn$Pn  ~  Pn$||HS(V*  ,H)  +  ||Pn $  ~  $ || HS( V* ,H) 

<  ||^Pn  ~  $|(HS(V,H)  +  |(Pn^  ~  ^||hS(V,H). 

Now  $  6  HS(V*,H)  implies  that  e  HS(H,V*)  and  that  ($P„)*  =  P^*  e  HS(H,V*) 
where  P*  denotes  the  adjoint  of  Pn  as  an  element  of  £(V*).  Then 
(4.5)  ||Pn^Pn  —  $||tfS(  V\H)  <  ||.P*$*  —  $*||//S(//,V)  +  ||Pn  $  ~  $||  HS(  V*  ,H)  • 

An  applicaiton  of  Lemma  4.2  with  X  =  H  and  Y  =  V*  and  the  strong  convergence  of  P* 
to  the  identity  on  V*  yield  that  the  first  term  on  the  right  hand  side  of  the  estimate  (4.5) 
above  tends  to  zero  as  n  — >  oo.  Similarly,  Lemma  4.2  with  X  =  V*  and  Y  =  H  and  the 
strong  convergence  of  Pn  to  the  identity  on  H  implies  that  the  second  term  tends  to  zero  as 
n  — ►  oo  as  well.  A  similar  argument  can  be  used  to  show  that  |jPn<t>Pn  —  $||hS(H,V)  ~ >  0 


as  n  — *  oo  and  the  lemma  is  proved. 


* 


I 

■ 


1 


Theorem  4.1  If  IIo  e  C(0,  T;  H)  is  the  unique  integral  solution  to  the  initial  value  problem 
(3.1),  (3.2)  (or,  equivalently,  (3.4),  (3.5))  and  II,,  e  C(0,T;Wn)  is  the  unique  integral 
solution  to  the  initial  value  problem  (4.3),  (4.4),  then  IIn  converges  to  IIo  in  C(0,  T;?f)  as 
n  — ►  oo.  That  is,  lim  |IIn(t)  —  IIo(t)| «  =  0  with  the  convergence  uniform  in  t  for  t  e  [0,  T], 

n  — ►oo 

Proof  The  desired  result  will  follow  immediately  from  Theorem  2.3  once  we  have  verified 
that  the  hypotheses  (i)-(iv)  given  in  the  statement  of  that  theorem  hold.  If  $o  e  Co  then 
<$„  =  Pn $0 P„  e  Cn  and  Lemma  4.3  implies  that  lim  |$n  —  $o|-H  =  0.  Thus  lim  Cn  D 

n  —►oo  n — oo 

Co-  Lemma  4.3  also  implies  that  lim  |Fn(t)  —  Fo(t)[x  =  lim  |PnF0(t)Pn  —  Fo|-h  =  0 

n— *oo  n — ►oo 

for  almost  every  t  t  (C,T).  Properties  of  Hilbert-Schmidt  operators  and  the  fact  that  P„  is 
an  orthogonal  projection  yield  |Fn(t)|-H  =  |PnFo(t)Pn|-H  <  |F0(t)|w  e  Lj(0, T),  for  a.e. 
t  e  (0,  T).  Consequently  hypothesis  (ii)  is  satisfied.  Once  again  from  Lemma  4.3  we  obtain 
lim  |II°  -IlJJlw  =  lim  |PnIIgPn  —  IIo]*  =  0.  The  verification  of  hypothesis  (iv)  is  all 

n — ^oo  n — >oo 

that  remains. 

Let  A  >  0  and  let  <&„  e  C„,  n  =  0,1,2,...  with  lim  |$n  —  ^o|-H  =  0.  Set  =  Jn(A)$n, 

n — ►oo 

n  =  0,1,2,...  where  Jn(A)  =  (I  +  AAn)-1.  Then  recalling  (3.3)  we  have 


-Pn'I'oPoilv  <  AlLoi^n-Pn^oPul^n-Pn^oPnl^ 

=  [(I  +  AAn)^n  -(1  + AA0)^oPn,'I'n  -Pn^oPn]^  “  [^n  “  ^oPn ,  “  Pn  ^oPn]« 

+A  [L0{#o  -  Pn^oPn},^  ~  Pn^oPnlx  “  A[Fn('&n)Pn  ~  *0  )Pn ,  "  Pn  *0  P„]« 

=  [$„  -  P„$0Pn,  “  Pn^oP..]*  -  [¥,,  ~  Pn^oPn,  ~  Pn^0Pn]H 

+A[L0{1'o-Pn^oPn}Pn,^n  -Pn^oPnlx  -A[F0(^„)Pn  -^(Pn'foPnJPn, 

*n  -  Pn'I’oPn]*  +  A  [F0(V  0  )Pn  ”  ^o(Pn*0Pn)Pn ,  ~  P„*0Pn]K 
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<  [$B  -P„$0Pn,^n  -Pn^oPn]W  +  A[L0{^0  ~  Pn $oPn  }Pn ,  -  Pn'I'oPn]* 

+A[^0(^o)Pn  -^o(P|,^oPn)Pn,^n  -Pn^oPJw 

<  ||*„  -  Pn^oPnHv*  Il'Pn  -  Pn*0P«||V  +  A  llLo{^0  -  P»*oP„  }  P„  ||  y  ' 

\\*n  -  Pn^oPn||V  +  A  ||^0(^o)P»  ~  Wn  *oP„  )Pn  ||  y.  ||*n  ~  P„*oPn||V  . 

where  in  the  above  estimate  we  have  used  the  facts  that  [^o('J'n)  —  •?ro(Pn^oPn)Pn> — 
Pn^oPn]  >  0  and  P„  nonnegative  self-adjoint  and  compact  (being  of  finite  rank)  imply 
[Toi^nWn  -^o(Pn4,oPn)Pn,^n  -Pn^oPn]«  >  0-  Thus 

||*„  -  PnM'oPnHv  ^  -t  l$»  "  Pn*oP„|* 

-f(l/o/){((^4*(|£(v  v»)  |Pn|c(H)  +  IMIic(V,V*)  l-Pnlc(V)}ll^  Pn^oPnllv 

+  Z  l^o(^o)  ~  ^b(Pn^oPn)|w 

for  some  positive  constant  k.  Now 

|*»  ~  Pn$0Pn\-H  <  l$»  "  $oj*  +  “  P*$oPJ* 

and 

ll'^*ll£(V,V*)  |P«l£(H)  +  ll^llc(V,V)  \Pn\ciV)  —  +  |P«lc(V)) 

which  is  uniformly  bounded  in  n.  The  assumption  that  lim  <£n  =  4?o,  Lemma  4.3,  and 

n— < *oo 

lim  |tf0  —  Pn^oPnlw  =  0  together  with  the  continuity  of  the  map  0  — >  .Fo(0)  from  H 

n  — *oo 

into  itself  yield  lim  ||’I'n  —  P„1i,pPll||y  =  0-  Consequently  the  estimate 

n  — »oo 

lim  |^n  -  'J'oltt  <  K  lim  (|^n  -  ’J'oJIv 

n— *-oo  n— *00 

<  I<  lim  ||^„  -  Pn^oPnilv +  K  lim  ||Pntf0Pn  -  tfollv 

n— *oo  n— *oo 

together  with  Lemma  4.3  yield  the  desired  result. 
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Remark  4.1  Although  we  are  unable  to  demonstrate  that  the  hypotheses  of  the  existence 
result,  Theorem  2.1,  are  satisfied  when  the  initial  value  problem  (3.4),  (3.5)  is  considered 
in  the  space  V,  we  can  show  that  hypotheses  (i)-(iv)  of  Theorem  2.3  are  in  fact  satisfied 
in  the  stronger  V  —  topology.  More  precisely  if  IIq  e  V  and  Fo  £  L](0,T;  V),  then  Lemma 
4.3  is  sufficient  to  obtain  the  convergence  of  IIn  to  IIo  in  C(0,T;  V).  In  the  case  of  a  linear 
dynamical  system  (i.e.  when  Ao  £  £(V,  V*))  it  is  in  fact  possible  to  apply  Theorem  2.1  in 
the  space  V  and  therefore  obtain  a  complete  approximation  theory  in  V  (see  [1]). 

5.  An  Example 

We  illustrate  the  application  of  our  approximation  results  with  a  linear-quadratic 
optimal  control  problem  (see  [10])  involving  a  one  dimensional  heat  equation  (exten¬ 
sion  to  higher  dimensions  is  straightforward).  Let  H  =  L2(0, 1)  endowed  "with  the  usual 

l 

inner  product, (^,  ip)  =  /  V^t?)^7?)^'/-  Consider  the  optimal  control  problem  of  finding 

o 

u  e  L2((0,T)  x  (0, 1))  which  minimizes  the  quadratic  performance  index 

T 

J(u)  =  (£x(T,  -),x(T,  •))  +  /  (Qx(t,-),x(t,-))  +  r(u(t, -),u(t,-))dt 

o 

subject  to  the  linear  distributed  dynamical  system 

(5.1)  |f(t,t?)  -  ^a(?7)|^(t,77)  =  bu(t,J7),  t  >  0,  0  <  rj  <  1, 

(5.2)  x(t,0)  =  0,  x(  t,  1 )  =  0,  t>0 

(5.3)  x(0, 7/)  =  x°(??),  0  <  rj  <  1. 

where  a  e  Loo(0,l)  with  a(?/)  >  a  >  0,  a.e.  i]  e  (0,  1),  b,r  e  R  with  r  >  0,  x°  e  L2(0, 1), 
Q,QeCo  =  {$  £  HS(L2(0,1),L2(0,1)),$  =  >  0}  Note  that  Q.  Q  (  C0  implies 
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that  {Q<p)(ri)  =  f  g(rj,t)V>({)d£  and  (Q^)(»?)  =  /  q(rj,€)<ptf)d(  with  g,q  e  L2{{ 0,  l)x(0,l)), 
o  o 

g(rht)  =  $(£.»?),  9(V,Z)>  0,  q{rt,S)  =  q{£,ri),  >  0,  a.e.  (r?,  £)  t  (0, 1)  x  (0, 1 ). 

If  we  set  V  =  Hg(0, 1)  endowed  with  the  usual  inner  product,  <  tp,y  >  = 

i 

J  D(p(r/)Dfl’(rj)dt>,  and  corresponding  induced  norm,  ||  •  ||,  then  V*  =  H_1(0, 1).  Define 
o 

the  operator  A  e  £(V,V*)  by  {Ap>){il')  =  (aDi^.DC),  p>,  </'  e  V.  Then  (Ay>,^>)  > 

a||<,?||2,  t  V  and  the  restriction  A  of  the  operator  A  to  the  set  Dom(A)  =  {p>  e  V  : 
Ap>  f  H)  (=  H2(0, 1)  fl  Hg(0, 1)  when  a  is  sufficiently  smooth)  is  positive,  self-adjoint  and 
the  infinitesimal  generator  of  a  uniformly  exponentially  stable  analytic  (parabolic)  semi¬ 
group  of  bounded,  self-adjoint  linear  operators  (T(t)  :  t  >  0}  on  H.  Let  U  =  L2(0, 1) 
and  define  the  operator  B  e  £(U,H)  by  (Z?v)(?;)  —  b v(r/ ),  v  e  U,0  <  <  1,  and  define 

71  c  C(U)  by  (7lv){r])  =  rv(?/)i  v  e  U,0  <  r;  <  1. 

The  solution  to  the  initial  value  problem  (5.1)  -  (5.3)  is  given  by 

t 

x(t,-)  =  T(t)x°  +  /T(t-s)5u(s)ds,  t  >  0. 

o 

where  for  u  e  L2((0,T)  x  (0,  1))  we  have  used  the  shorthand  notation  t  — ►  u(t)  to  denote 
the  function  t  — +  u(t,  •)  e  L2(0,T;  U).  The  solution  to  the  optimal  control  problem  is  given 
in  closed- loop,  linear  state  feedback  form  (see  [10])  by 

(5.4)  u(t)  =  —  'R-1£>*n0(T  —  t)x(t,  •)  =  —  (b/r)n0(T  -  t)x(t.  •),  a.e.  te(0,T), 

where  n0  is  the  nonnegative  self-adjoint  solution  to  the  Riccati  differential  equation  initial 
value  problem 

(5.^)  n0(<)  +  >tn0(t)  +  n0(tM  +  rio(t)iJ7v  3 ^ * n 0 ( t )  =  Q  a.e.  t  e  (0,  T) 
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(5.6) 


IIo(O)  =  Q. 


If  we  set  M  =  BTl~lB*  =  (b2/r)I  e  £(H),  then  M  =  M*,M  >  0  and,  since  $M  = 
e  C0,  >  0,  $  e  C'0.  Thus  Temam’s  [15]  theory  (with  F0(t)  =  Q,  t  e  [0,T], 

=  (b2/r)<$2,  and  Ily  =  Q)  presented  in  section  3  above  yields  the  existence 

of  a  unique  integral  solution  ITo  e  C(0,  T;  HS(H,  H ) )  to  the  initial  value  problem  (5.5), 
(5.6)  with  n0  e  L2(0,T;  HS(H,  V)  n  HS(V*, H)),  ri0  e  L2(0,  T;  HS(H,  V*)  +  HS(V,H)) 


and  n0(t)  e  Co,  t  e  [O.T].  (If  Q  is  such  that  AQ  +  QA  e  HS(H,H)  -  for  example 
if  Q  =  0  -  then  IIo  will  be  a  strong  solution  with  IIo  e  W1,oo(0,  T;  HS(H,  H))).  Since 
H  =  HS(L2(0, 1)),  L2(0, 1))  is  isometrically  isomorphic  to  L2((0, 1)  x  (0, 1)),  it  follows  that 
there  exists  a  7r0  e  C(0,T;  L2((0, 1)  x  (0, 1)))  with  7r0(t,»7,^)  =  tt0 ( t ;  £ ,  r? ) ,  7r0(t,77,f)  > 
;  0.  a.e.  It]  £)  e  (0, 1)  x  (0. 1 ) ,  t  e  [0,T]  such  that  the  solution  to  the  optimal  control  problem, 

(5.4),  is  given  by 

i 

u(t,?/)  =  --  J  tt0(T  -  t;  ?7, £)x(t, £)d£, 

I 


for  almost  every  (t,^)  e  (0,T)  x  (0,1). 

We  consider  a  linear  spline  based  approxirnaton  scheme.  For  each  n  =  2,3,...  let 
Hn  =  span{^i,  Jj'Tj1  where  for  j  =  l,2,...,n-l,  i^n  is  he  jth  linear  spline  ("hat”)  function 
on  [0,1]  defined  with  respect  to  the  uniform  mesh  {0, 1/n,  2/n, ...,  1}.  That  is 


VJn(V)  =  { 


0  0  < q  < 

nn-j  +  1  ^r<V<i; 

j  +  f  i  <  n  <  2^i 

0  *—<»/<!• 


Let  Pn  :  H  — >  Hn  denote  the  orthogonal  projection  of  H  onto  H„  with  respect  to  the 
(•,•)  inner  product  and  define  ,4,,  e  £(H„)  to  be  the  Galerkin  approximation  to  A.  More 
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precisely,  we  set  —  Vn,  e  Hn,  where  yn  is  the  unique  element  in  Hn  which 

satisfies  (A?n,9„)  =  {i'a,9n),  9n  e  Hn.  Let  Qn  =  P„Q  e  £(Hn)  and  Qn  =  Pn(?e£(H„). 

Using  the  properties  of  interpolatory  splines  and  density  arguments  it  is  not  difficult  to 
show  that  lim  ||P„<^>  —  v? I j  =  0,  e  Hg(0,l)  =  V  and  consequently  that  the  assumption 

n— oo 

(4.1)  is  satisfied.  It  follows  therefore,  that  the  approximation  theory  developed  in  section 
4  yields  that  the  solution  IIn  to  the  finite  dimensional  Riccati  differential  equation  initial 
value  problem 

(5.7)  rin(0  +  An„(t)  4-  nn(t)A  +  (b2/r)nn(t)2  =  Qn,  a.e.  t  €  (0.T), 

(5.8)  nn(0)  = 

with  nn(t)  nonnegative  and  self-adjoint,  satisfies  lim  |IIn(t)Pn  —  XTo ( t ) I HSf H  H)  =  0 
uniformly  in  t  for  t  e  [0,  T] . 

Since  the  basis  {^n}j*T11  is  not  orthonormal,  simply  replacing  the  operators  in  (5.7). 

(5.S)  with  their  matrix  representations  will  not  lead  to  the  familiar  symmetric  matrix 

Riccati  differential  equation.  For  a  linear  operator  Ln  with  domain  and/or  range  in  Hn. 

we  denote  its  matrix  representation  with  respect  to  the  basis  by  L\.  Define 

:  [0. 1]  -*•  Rn_1  by  $1,(77)  =  (v’nC^), ---- (^ ))T  and  set  ,WN  =  (<$„•$„)  = 
1 

J  $n(ri)$n(n)Jdi].  Then  Av  =  (aD$n,  ),  G.y  =  $„)  and  Qn  = 

0 

If  we  bt  fe  =  M n Gs  1  Qn  =  M'sQ^i  and  II x ( t )  =  XNnx(t).  then 
nN  is  the  solution  to  the  initial  value  problem 

* 

n.v(t)  +  -4X fls ( t )  +  n\(t)«4\  +  (b2/r)  nx(  t  ).Vf  ^ 1 17x  ( t )  —  Qs •  t  f  ( 0.  T ) 

n.v(Q)  =  Gx 


25 


with  ITn (t )  nonnegative  and  self-adjoint.  The  approximating  solution  to  the  optimal  con¬ 
trol  problem  then  takes  the  form 

1 

un(t,r?)  =  "  j  ?rn(T  -  t;r/,£)x(t,£)d£,  a.e.  (t,  ?/)  e  (0,  T)  x  (0, 1) 
o 

where  7rn  e  C(0,T;  L2((0, 1)  x  (0, 1)))  is  given  by 

7Tn(t;T/,0  =  $n(>?)T^N1nN(t)A<^1$n(^)1 

for  t  e  [0,  T],  (?;,  £)  e  [0, 1]  x  [0, 1].  Our  convergence  theory  yields  lim  7ru  =  7r0  in  L2((0, 1)  x 

n— *oo 

(0, 1))  uniformly  in  t  for  t  e  [0,T].  That  is 

l  i 

lim  sup  /  /  |7r„(t;r?,£)  -  7r0(f;  r/,  £)|2  dj]d(,  =  0. 

n-~°°  <{[o,T1  J  J 
0  0 

We  note  that  some  interesting  open  questions  which  are  related  to  our  effort  here 
remain.  For  example,  a  Hilbert-Schmidt  existence,  approximation,  and  convergence  theory 
for  the  case  in  which  the  operator  A  (and  therefore  the  operator  A0)  is  time  varying.  The 
standard  existence  and  approximation  results  in  the  spirit  of  Theorem  2. 1-2.3  for  nonlinear 
systems  involving  time  dependent  operators  are  typically  too  restrictive  to  be  of  any  use 
in  the  Hilbert  space  framework  in  which  the  problem  would  be  formulated.  An  approach 
such  as  the  one  which  was  taken  in  [2]  in  the  context  of  approximation  methods  for  inverse 
problems  for  nonautonomous  nonlinear  distributed  systems  which  was  based  upon  Barbu's 
[3]  Theorem  III.  4.2  may  be  more  appropriate. 

Finally  a  Hilbert-Schmidt  approximation  theory  for  the  steady  state,  or  algebraic 
operator  Riccati  equation  can  also  be  developed  (see  [12]).  (The  algebraic  Riccati  equation 
arises  in  the  context  of  the  linear-quadratic  optimal  control  problem  on  the  infinite  time 


interval). 
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